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Abstract. In this paper continuity theorems are estabhshed for the num- 
ber of losses during a busy period of the M/M/l/n queue. We consider an 
M /GI /1/n queueing system where the service time probabiUty distribution, 
slightly different in a certain sense from the exponential distribution, is approx- 
imated by that exponential distribution. Continuity theorems are obtained in 
the form of one or two-sided stochastic inequalities. The paper shows how 
the bounds of these inequalities are changed if further assumptions, associated 
with specific properties of the service time distribution (precisely described in 
the paper), are made. Specifically, some parametric families of service time 
distributions are discussed, and the paper establishes uniform estimates (given 
for all possible values of the parameter) and local estimates (where the param- 
eter is fixed and takes only the given value). The analysis of the paper is 
based on the level crossing approach and some characterization properties of 
the exponential distribution. 



1. Introduction 

1.1. Motivation, the problem formulation and review of the related lit- 
erature. In a large number of engineering applications of queueing theory it is 
pertinent to know what to expect when we replace one probability distribution 
function of a given stochastic model, by another probability distribution with sim- 
pler/concise properties. Many of these applications are associated with the case 
where a probability distribution is replaced by an exponential distribution if it is 
known that the given distribution is close (in a certain sense) to the exponential 
distribution. These problems are closely related to characterization problems asso- 
ciated with the exponential distribution and its stability (e.g. Azlarov and Volodin 

Continuity analysis of queueing systems is a difficult problem when compared 
to general characterization and continuity problems for random variables, and it is 
relevant to know the expected behavior of the system under a range of different con- 
ditions. The paper discusses three different conditions under which the probability 
distribution function of a service time, which differs slightly from the exponen- 
tial distribution in the uniform metric, can be approximated by that exponential 
distribution. 

This paper considers an M /GI /1/n queueing system, where n is the buffer size 
excluding a customer in service, A is the rate of Poisson input, B{x) — Pr{x < x} 
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is the probability distribution of a service time Xi a-nd the parameter fi is the 
reciprocal of the expected service time. It is assumed then that the probability 
distribution function B{x) is close (under specific cases described below) to the 
exponential distribution with the same parameter /z. 

Let Ln denote the number of losses during a busy period of this queueing system. 
The classic explicit results for the M/GI/l/n queueing system, for example the 
recurrence relations for the expectations EL„, as well as a number of results on the 
asymptotic behavior of EL„ as n — > cxd, can be found in Abramov [I], [1], Cooper 
and Tilt [T7^ and Tomko [35 . An application of these explicit results to asymptotic 
analysis of lost messages in communication networks is given in Abramov |9j , and 
their application to problems of optimal control of large dams is given in [llj and 
[12j . Some simple stochastic inequalities for different finite buffer queueing systems 
are obtained in Abramov [5], [TOj, Righter [32], Pekoz et al. [29], Wolff [^ and 
other papers. Other relevant studies of the loss probability in finite buffer systems 
can be found in Abramov , f8] and Choi et al. [M] . 

The present paper establishes continuity theorems in the form of stochastic in- 
equalities for Ln similar to the stochastic inequalities established in Abramov [5], 
[lOj . The main difference between the results of this paper and those of Abramov 
0, [in] is that the stochastic inequalities established in the present paper depend 
on specific properties of the service time distribution, whereas the stochastic bounds 
of the earlier papers [5] , |10| are obtained without any special assumption on that 
service time distribution. Then the specific feature of the present paper is that 
the stochastic inequalities obtained can be considered as approximation bounds or 
continuity bounds for the queueing system in the uniform metric, and these approx- 
imation bounds can have real application for stochastic systems arising in practice. 
The phrase continuity theorems used in this paper is associated with the traditional 
" "e-(5 language" of mathematical analysis" . The input characteristic of the system 
is assumed to have a small variation e, and then the output characteristic has a 
small variation 5 continuously depending on e. Both these e and 5 variations are 
assumed to be described by some special metrics. A large variety of metrics for the 
continuity analysis of complex systems can be found in Dudley J.9J, Kalashnikov 
and Rachev [26], Rachev [30], Zolotarev [40] and many others. In the cases where 
an application of the well-known traditional metrics such as the uniform metric or 
the Levy metric becomes unavailable, very difficult or not profitable, the continuity 
analysis of a complex system requires to use the special type of metrics, appropriate 
for analysis of a given system or associated with given special conditions. 

The continuity analysis of queueing systems goes back to the papers of Kennedy 
[27] and Whitt [36] and further developed by Kalashnikov [25], Zolotarev [38], [39] 
and Gordienko and Ruiz de Chavez [22], [23] to mention a few. 

The continuity analysis of such standard characteristics as the number of losses 
during a busy period is a very difficult problem. One of the advantages of the 
method in this paper, is that the continuity of this characteristic is studied in 
terms of the uniform probability metric. As a result, the final results are simple 
and clear, and the conditions of the theorems are verifiable. (The uniform metric 
of one-dimensional distributions is usually called the Kolmogorov metric.) Further- 
more, under different assumptions that the service time distribution belongs to its 
parametric family of distributions depending on parameter p we establish uniform 
estimates for the number of losses during a busy period given for all possible values 
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p. That is, stochastic inequahties depending only on e but not on p are obtained. 
In the case where parameter p of this family is fixed (specifically the case P ~ \ 
is considered) we establish essentially stronger estimates than these corresponding 
uniform estimates. 

For the purpose of the continuity analysis, this paper essentially develops the 
level crossing approach of the earlier papers of Abramov [T, f3], f5], f6], [TCF and 
some characterization theorems associated with the exponential distribution, that 
are then used together with the aforementioned metrical approach. We also develop 
the continuity theorem of Azlarov and Volodin fl3' under some additional conditions 
and use the properties of the known classes of special probability distributions such 
as NBU (New Better than Used), NWU (New Worse than Used) (e.g. Stoyan 
[51] for the properties of these classes of distribution and related results). The 
special parametric order relation which was introduced and originally studied 
in Abramov [[IJ, is also used. For the reader's convenience all necessary facts and 
concepts are recalled in this paper. 

1.2. The main conditions for the probability distribution of a service 
time. The main conditions for the probability distribution function B{x) (i?(0) — 
0) of a service time under which we obtain the stochastic bounds are as follows. 

• Condition (A). The probability distribution of a service time has the represen- 
tation 

(1.1) B{x) = pF{x) + (1 - p)(l - e-''^), < p < 1, 

where F{x) — PrK' < x} is a probability distribution function of a nonnegative 
random variable having the expectation ^ , and 

(1.2) sup \Fy{x) ~ F{x)\ < e, e > 0, 

a:,y >0 

where Fy{x) = Vr{C, < a; + y | C > v}- Relation (II. 2p says that the distance in 
the uniform metric between F{x) and E^{x) = 1 — c^^^ is not greater than e, the 
case e = corresponds to the equality F{x) ~ Efj_{x) for all x. In the following 
we assume that e > and write the strong inequalities (i.e. the right-hand side of 
(|1.2p is less than e). 

• Condition (B). Along with (|l.ip and (|1.2p it is assumed that F{x) belongs 
either to the class NBU or to the class NWU. 

Recall that a probability distribution function S(a:) of a nonnegative random 
variable is said to belong to the class NBU if for all a; > 0, y > we have S(a; + 
y) < S(a;)S(2/). (Throughout the paper, for any probability distribution function 
S(x) we use the notation S(a:) = 1 — 5(a;).) If the opposite inequality holds i.e. 
S(a; + y) > 'E.{x)E{y) then S(a;) is said to belong to the class NWU. 

• Condition (C). The probability distribution function Bix) of a service time 
belongs to the class NBU, and 

(1.3) sup \By{x) - B{x)\ < e, 

3:,i/>0 

where By{x) — Pr{x < x + y \ x > u}- 

• Condition (D). Let Fy{x) = Pr{( <x + y\ (>y}hea family of given 
probability distributions. It is assumed that there exists the probability distribution 
Fyo{x) of this family satisfying the relation FyO <<^^ Fy for all y >0. 
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The definition and main property of tlie parametric order relation "^x is recalled 
in Appendix A. For a more detailed consideration see Abramov [1]. 

• Condition (E). Let By{x) — Pr{x <2; + ?/|x>2/}bea family of given 
probability distributions. It is assumed that there exists the probability distribution 
Byo{x) of this family satisfying the relation ByO By for all y > 0. 

If the family of probability distributions Fy{x) (or By{x)) is partially ordered 
with respect to the order relation '^\, then according to Zorn's lemma (e.g. Ciesiel- 
ski [15^ ) FyQ (or correspondingly ByO ) is a minimal element of this family (regarding 
this order relation 'lox). 

The only difference between Conditions (D) and (E) is that Condition (E) is 
related to the probability distribution function B{x) of a service time, whereas 
Condition (D) is related to the associated probability distribution function F{x). 

Conditions (A), (B) and (C) are the main conditions for our consideration, while 
Conditions (D) and (E) are additional (associated) conditions. In other words one 
of Conditions (A), (B) and (C) are always present in our consideration. Condition 
(D) can be present only with Conditions (A) and (B), and Condition (E) can be 
present only with Condition (C). 

It is worth noting that in certain known cases the minimal element FyO can 
be determined easily. For example, in the case where F{x) belongs to the class 
NWU, the minimal element of the family of probability distributions Fy{x) is FyO = 
Fq — F va the sense of the order relation . (By replacing NWU with NBU and 
respectively 'minimal' with 'maximal', this fact is explicitly used in Theorem 14.21 
Section 4.) 

1.3. Further discussion of the main conditions. Conditions (A), (B) and (C) 
are three different conditions where if e is small, then the probability distribution 
function B{x) is close to the exponential distribution in the sense of the uniform 
metric. Therefore it is very significant to know what one can expect if Conditions 

(A) , (B) or (C) are satisfied. 

In most queueing problems representation (|1.1|) is not standard. The standard 
assumption appearing in characterization problems associated with the exponential 
distribution is (jl.Sp (e.g. Azlarov and Volodin [T3j; cf. Daley [18 for some earlier 
results). Representation of the probability distribution B{x) in its special form 
(|l.ip matches the usual case B{x) = F{x) (that is the case p = 1). It is shown 
below (see Lemma [2T|) that assumptions (|1.1|) and (|1.2p of Condition (A) allow us 
only to prove that for all < p < 1 

(1.4) sup \By{x)-B{x)\ < 5e, 

3:,y>0 

and only in the special case where p = ^ we have 

(1.5) sup \By{x) -B{x)\ < 2e 

(see Remark 12. 3p . As a result, the approximation bounds under Conditions (A), 

(B) and (C) are all different, and more specifically representation p.ip together 
with assumption (|1.2[) lead to relatively worse estimates than those in the special 
case B{x) = F{x) under assumption (|1.3|) . However (|1.4p is a uniform estimate 
for all < p < 1, whereas (jl.Sp is a usual local estimate better than (|1.4p . In the 
special case p = 1 it is a local estimate too, which in fact is (|1.3p . 
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The advantage of is also as follows. When the parameter p is small, the 
service time distribution is close to the exponential one in the uniform metric. 
Thus, our results are associated with a wider class of approximations and based 
on a two-parameter family of distributions. In practice, this can help us find a 
more appropriate representation and approximation for the initial probability dis- 
tribution function B[x) given empirically. The main results of this paper establish 
explicit dependence of the parameter e only and are uniform in p. However there 
are examples related to the special case oi p — ^ where the bounds obtained are 
essentially better than in the general case of arbitrary p. In a similar manner the 
case of any given p can be considered, where the specific estimates are expected to 
be better than the corresponding uniform estimates. 

Notice also that the class of stochastic inequalities is wider. On the one hand, 
only the parameter e can be assumed to be small, and on the other we can assume 
that both e and p are small values. The variety of these assumptions enables us to 
choose the possible appropriate parameters of the model for the further approxi- 
mations in order to obtain a relevant conclusion. 

1.4. Organization and methodology of this paper. This paper is structured 
into 4 sections. Following this introduction, Sections 2, 3 and 4 study the number 
of losses during a busy period under Conditions (A), (B) and (C), respectively. 

In Section 2 we give a deeper analysis of the intervals obtained under the special 
procedure of deleting subintervals and connecting the ends as is explained in a 
number of the earlier papers of the author (see Abramov [1], [3], [5], [6], [10]). 
Specifically we prove the following results. The first result of this section, Lemma 
12.11 provides the uniform estimate (|1.4p for the parametric family of distributions. 
Remark 12.31 provides local estimate (|1.5|) in the special case p = ^. The proof 
of both these estimates is based on the characterization theorem of Azlarov and 
Volodin [13j, which is formulated in Lemma [2T2l The proof of Theorem 12. 41 is based 
on the level crossing approach which in its present form is originated and developed 
by the author (see Abramov [I]. [3]. [5]. [6]. [lO]). 

In Section 3, under Conditions (B) and (D) stronger results than those under 
Conditions (A) and (D) are established. The results are based on similar proofs to 
those in Section 2. 

In Section 4, under Condition (C) the two-sides stochastic inequalities are given 
in Theorem 14.21 The proofs of the statements under Condition (C) are completely 
analogous to the earlier proofs under Condition (A). 

The partial order relation and its main property are discussed in Appendix 
A. The proofs of Lemma 12.21 and Lemma |3. II is given in Appendix B. In Appendix 
C, the proof of Lemma [2?6l is given. 

2. The number of losses under Condition (A) 

Throughout the paper we use Kolmogorov's (uniform) metric between two one- 
dimensional probability distributions. Recall the definition of Kolmogorov's metric 
(e.g. Kalashnikov and Rachev [55], Rachev [3D]). Let G{x) — Pr{^ < x} and 
H{x) = Pr{?7 < a;} be probability distribution functions of the random variables ^ 
and rj. Kolmogorov's metric J^{G, H) between two probability distribution func- 
tions G{x) and H{x) is defined by 

JT (G, H) ^ sup I G{x) - H{x) \ . 
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In the following for Kolmogorov's metric the notation J(^{£,, if) is also used, where 
■^{^tV) — ■^{G,H). All exponential distributions are denoted Ea — Ea{x) — 
1 — e~"^, where a > is the parameter of distribution. 

Lemma 2.1. Under assumptions and il.'2^) we have 

sup \By{x) — B{x)\ < be. 

x>0,y>0 

In order to prove this lemma we use the following result of Azlarov and Volodin 
[13j . the proof of which is also provided in Appendix B. 

Lemma 2.2. (Azlarov and Volodin [13].^ Let F{x) = Pr{C < x} be a probability 
distribution function of a nonnegative random variable with the expectation 
and 



sup 

x>0,y>0 



\Fyix) - Fix)\ < e. 



Then 



je{F,E^) < 2e. 

Proof of Lemma 2.1. We have 

I - pF{x + y) - [l - p)E^{x + y) 



Byix) 



(2.1) 



Therefore, 



l-pF{y)-{l-p)E,,{y) 
pF{x + y) + (l- p)e-i'i'''+v) 
pF{y) + {1 - p)e-t'y 



(2.2) 



\Byix)-B{x)\ 

pF{x + y) + {l- p)e~^(^+2') 



pF{y) + (1 - p)e-f y 
pF{x + y) + (1 - p)e-''(^+2') 



-pF{x) - (1 -p)e-''" 
P^F{x)F{y) 



pF{y) + (1 - p)e-f y pF{y) + (1 - p)e^^'y 

p{l - p)[e-^''■'F{y) + e-^'yF{x)] (i ^ pfe-t^i^+y) 



(2.3) 



pF{y) + (1 -p)e-^^y 
According to (|1.2p we have 

F{x + y) 



F{y) 



F{x) 



pF{y) + il-p)e- 



where the denominator in the fraction of (|2.3p is assumed to be positive, and if it 
is equal to zero, then we use the convention ^ = 0. It follows from Lemma that 



(2.4) 



lF(a;) -e-^^1 < 2e. 
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Therefore from (l2?2l) . ([Qj) . (12^ and the triangle inequality for all a: > 0, ?/ > 0, 
we obtain: 

\pF{x + y)- p'F{x)F{y) - p(l - p)e-'"F(y) I 



< VF{V) 



: pF{y) 



F{x + y) 



-fix 



_F{y) 

F{x + y) 



+ p^F[y) F{x) -e 



fix 



(2.5) 



< 



-p''F{y)\F{x) 



F{x + y) 



F{x) + F{x) 



-fix 



+ p^F{y)\F{x)- 



- F{x) + \ F{x) -e" 

-fix I 



fix 



< iepF{y) + 2ep'F{y) 

< 5epF(y), 

I (1 - p)e-^(^+^) - p(l - p)e-''yF{x) - (1 - p)2e"^(^+^) I 



e^^^ - F{x) 1 + 1(1- p)F{x) - (1 - p)e-^^ | 
2e(l -p)2e-^^ 



< (1 -]5)e-^^ 
(2-6) <2e(l-p)c-^^ 

< 3e(l - p)e-^'y + 2e(l - pfe-^'y 

< 5e{l-p)e-''y, 

and then it follows from p.5p and p.6p that for all x > 0, y > 0, 

ii.,(.)--i.(.)i< ^-^g;^;+^-(^7)-""^ ^5.. 

Hence, Lemma |2. II is proved. 

Remark 2.3. The statement of Lemma |2. II is a uniform estimate for all p. In special 
cases where p is given, one can obtain stronger estimates. For example, in the case 
p = i we have the following inequalities: 



|iF(x + y) ^ lF{x)F{y) - le-'^-Fiy)\ 



(2.7) 



< jFiy) 


F{x+v) 
F(y) 


F{x)\ 


+ 4^^y)\ F{y) 




< \F{y) 


F(x+y) 

F(y) 


Fix)\ 


+ \F{y)\F{x)- 





< ieF(y) + ieF(y) 
- ^F{y), 



-fi{x+y) 



Fix) 



-fj.{x+y) I 



F{x)\ 



(2.8) 



< iee-'^'' 



< ee 
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Therefore, from (1271) and (l2?8l) we obtain 



(2.9) \By{x) - B{x)\ < ^ = 2e. 

Theorem 2.4. Under Conditions (A) and (D) we have 

Ln ^st 

1=1 

Z„ denotes the number of offspring in the nth generation of the Galton-Watson 
branching proeess with Zq ^ I and the offspring generating function 

, , B(\) l-z + 2zB{X) , , 

2.10) g{z) < + 5e — ^A^. z < 1, 

l-z + zB{X) l-z + zB{X) 

where 

B{X)=pF{X) + {l-p)^, 

fi ^ A 

F{X) = / c-^^di^(x). 
Jo 

I = 1, 2, . . is the sequence of nonnegative integer-valued independent identically 
distributed random variables all having the probability law 

(2.11) Pr{Q=m}=/ e-^^^^dG(x), 

Jo 

and G{x) is a probability distribution that satisfies the inequalities G <<g'^ FyO and 
G '^.'^x E^fj.! where FyO Fy for all y > 0. 

Remark 2.5. The condition (|2.1ip of Theorem l2.4l involvcs a probabihty distribution 
function G{x). If Condition (D) is not satisfied, or the probabihty distribution 
Fyo{x) is unknown and cannot be evaluated, then a probabihty distribution function 
G{x) should be chosen such that G fo'" all y > and G ^tJ.- 



Proof. Consider a busy period of an M/GI/l/n queueing system. Let f{j),0 < j < 
n+1, denote the number of customers arriving during a busy period who, upon their 
arrival, meet j customers in the system. It is clear that /(O) — 1 with probability 
1. Let ij,i,ij,2, tj.f(j) be the instants of arrival of these /(j) customers, and let 
Sj,!, Sj,2, Sj,f{j) be the instants of service completions (the case j < n) or losses 
(the case j = n + 1) at which there remain only j customers in the system. Note 
that tn+i^k = Sn+i,k, 1 < < f{n + 1), and f{n + 1) = L„, the number of losses 
during a busy period. 

For < j < n consider the intervals: 

(2.12) [tj,i,Sj^i), [ij-,2,Sj,2),---, 
It is clear that the intervals 

(2.13) [tj+i^i,Sj+i^i), [tj+1^2, Sj+1,2), 

are contained in intervals (|2.12p . Let us delete the intervals in (|2.13|) from those 
in (|2.12p and connect the ends, that is every point ij+i,fc with the corresponding 
point Sj+i,fc, 1 < fc < f{j + 1), if the set of intervals (|2.13p is not empty. In other 
words, in the interval of the form [tj,k,Sj^k), the inserted points are of the form 
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tj+i,m- The random variable ^j^fc is then the length of interval [tj^k, ■Sj.fc) niinus the 
intervals [ij+i,m, contained in this interval. 

A typical example of the level crossings on a busy period is given in Figure 1. 
The arc braces in the figure indicate the places of connection of the points after 
deleting the intervals. 

In the example, given in Figure 1, n = 3 and the number of inserted points of 
level 1 is equal to 2, that is /(I) = 2. Similarly, /(2) = 1. 

Let us take one of the intervals, [tj.k, Sj.fc) say, and a customer in service at time 
tj^k- Let Tj^k denote the time elapsed from the moment of a service start for that 
customer until time tj^k, and let i}j^k denote the residual service time of the tagged 
customer. The analysis of a residual service time is well-known in the literature, 
and in the simplest cases is associated with the inspection paradox or waiting time 
paradox (see e.g. [16], [20], [21], [23, [28], [31], [33] and many others). 

The sum r^jt + -dj^k has a more complicated structure than that in the stan- 
dard inspection or waiting time paradox, and deriving an explicit representation 
for Pi'l'dj k < x} in terms of the probability distributions B{x) or By{x) is very dif- 
ficult. By using Lemma 1 2 . 1 1 above . the probability distribution function Bj_k{x) — 
Pr{i?j_fc < x} can be evaluated nevertheless. 

Note first the following technical result. 

Lemma 2.6. For any integer j >1 and k >1, and any nonnegative real x and y, 

Pr{i?j,fc < X I Tj,k = y} = By{x). 
The proof of this lemma is given in Appendix C. 

Therefore, from Lemma 12.61 and Lemma l2. II for any nonnegative y we have: 

(2.14) \Pr{dj^k < X I T,- fc ^y}- B{x)\ = \By{x) - B{x)\ < 5e. 

We will show below, that by using the formula for the total probability one arrive 
at 

(2.15) |Pr{i?j,fc < a;} - B{x)\ < 5e. 

For this purpose we should show that the random variable Tj.fc is proper, i.e. 
PrjTj^fc < oo} — 1, and then use the formula for the total probability: 

/•oo 

Pr{i?,, fc < .t} = / Pr{i9j,fc < X \ Tj,k = 2/}dPr{T,- fc < y } 

(2.16) ^° 

= / By{x)dPi{Tj^k < y } 
Jo 

(Notice, that the random variables 'dj,k and Tj^k are generally dependent, and we 
cannot derive a representation for Pr{rj_fe ^ ?/}■) 

Let us show that the random variable Tj^k is proper for any j and k. Indeed, 
let Ai denote the event that there are i arrivals during the time interval from the 
service beginning of the customer, who is currently being served at the moment 
tj^k, until the time moment before the arrival at tj^k (the arrival at time tj^k and 



10 



VYACHESLAV M. ABRAMOV 




k,\ fl,1 h,\ S2.1 3|.| t\2 31.2 Sol 



Figure 1 . Up- and down-crossings on a busy period 
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possible arrivals at the moment of the service began are excluded). Then, 

(2.17) Pr{r,-fe <a;} = ^ Pr{T,, fe < x | A}Pr{A} 

i=0 

(2.18) > Prjsum of i + 1 interarrival times < a;}Pr{Ai} 

1=0 

(2.19) > Pr{sum of j interarrival times <x}. 

(The expression for the probability distribution of (j2.19p can be written explicitly.) 
Inequality (|2.18p follows from (|2.17p due to the following property: 

Pr{rj, fc <x\Ai} 

:= Pr{sum of i + 1 interarrival times < a; | all of them occur during a service time} 

> Pr{sum of i + 1 interarrival times < x}. 

Inequality (I2.19P follows from (|2.18p due to the elementary fact that for i < j 

Pr{sum of i interarrival times < x} 

> Pr{sum of j interarrival times < a;}. 

Thus, Tj^fe is a proper random variable, and (|2.16p is valid. Therefore, denoting 
the probability of (|2.16p by Bj^k{x), we see that under Conditions (A) and (D): 

\Bj^k{x) - B{x)\ < / sup \By{x) - B(x)| dPr{Tj- < y } 

Jo x>0,y>0 

(2-2'^) <5e due to Lemma [^TT] 

< 5e / dPr{Tj- fc < y } = 5e, 

and (|2.15p follows. 

Let Kj^k denote the number of inserted points within the interval [ij./t, Sj,fc) 
iEkll ^j.k := fU + 1)). Then for 1 < j < n (n > 1) 

Pr{K,,k^O}= / e-^M%fc(a;) =B,,fc(A), 
(2.21) Jo 

Pr{Ac,- fc = m} = [1 - Pr{^j,k = 0}][1 - B{X)r-^B{X), m = 1, 2, ... 

where B{X) ~ e~'^^dB{x). Notice, that according to (|2.20p for any 1 < j < n 
{n > 0) 



|B,,,(A)-B(A)| 



(2.22) 



/ e-^^dBj,fc(x) - / e-^'^dBix) 

Jo Jo 

< Xe-^^'lB.jix) - B{x)\dx <5e. 

Jo '■ V ' 



<5c due to iT20i 



The term 1 — Pr{Kj k — 0} of (|2.2ip is the probability that during the residual 
service time of the tagged customer there is at least one arrival (and therefore at 
least one inserted point). The next term [1 — B{X)]"^^^B{X) means that along with 
the first inserted point associated with the tagged customer, there are m — 1 other 
inserted points. 
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For j — n we have 

(2.23) Pr{«;„.fc = m} = / e-^^^-^dS„,fe(a;), m = 0, 1, ... 

Let us now consider the time interval [to.i? Sq.i) and the number of inserted points 
Ko,i- Taking into account that Bq^i(x) = Bix) we have 

(2.24) Pr{Ko,i = to} = [1 - ■§(^)]''"-B(A), to = 0, 1, 

in the case when the loss queueing system is not M/GI/l/O, and 

-\x (Aa;)™ 



-dB{x), m = 0, 1, 



in the case of the M/GI /I/O queueing system. 

In the case when the queueing system is not M/GI/1/0, for the probability 
generating function of Ko,i from (|2.24p we have the representation 

(2.25) y z™Pr{Ko 1 - m} = , |z| < 1. 

^0 ' l-z + zB{X) 

Relation (|2.25p has been established in [TD] as a particular case of a more general 
result. 

Let us now prove the inequality (|2.10p . For all j < n and k > 1 denote: 



9j 



fc(z) = ^ z™Pr{K,, fc = to}, |z|<L 



Comparing this expression with (|2.25p we have as follows. 
Let us denote 

B{X) 



Then from (P?^ and (P?^ we have 

\9*{z) - gjM{z)\ 



1- z + zB{X) 
^ z™ [Pr{Ko,i = m} - Pr{Kj,fc = to}] 

m=0 

oc 

B{X) ~ B,,fc(A) + ^ - B(A)]"-iB(A)[%,(A) - B{X)] 



< 



\B{X) - %fe(A)| + |B(A) - S,,fe(A)| z y z^[l - B{X)rB{X) 
V — ■ ' " V — ■ ' - 

<5e due to l lOa <5e duC tO iT22\ 

<5e[l + zg*iz)]. 



m=0 



Therefore, 



gj^kiz) <5e[l + zg*{z)]+g*{z) 



1- z + 2zB(X) 

^ + 5e 

l-z + zB{X) l-z + zB{X) 

for all j < n and fc > 1 , and we arrive at 

, s /X B(X) 1- z + 2zB(X) 

(2.26) g(z)= sup g,Az) < + 5^ • 

j<n,k>i 1 - z + zi3(A) 1 - z + zi3(A) 



CONTINUITY THEOREMS FOR THE M/M/l/n QUEUEING SYSTEM 13 

and (|2.10|) follows. Notice, that from the detailed calculations of p.6p . inequahty 
^(TM) is strong. 

Now, in order to finish the proof of the theorem, let us first establish some 
elementary properties of the class of probability distribution {By{x)}. Recall that 
the explicit representation for By{x) follows from (|2.ip . 



Property 2.7. The probability distribution By{x) can be represented as 

(2.27) By{x) = TyFyix) + (l ^ T y) E ^[x) , 

with parameter ry < 1 depending on y. 
The equivalent form of (|2.27p is 

By{x)^ryFy{x) + {l~ry)e-'^\ 

Proof. Indeed, we have: 

(F{x + y) ^ ^ „ _ pF{x + 2/) + (1 - p)e-^^i-+y) 



After some algebra we see that 

^ pF{y) 
^ pF{y) + {\-p)e-i^y' 
which clearly constitutes that ry < 1. □ 

We have the following property. 
Property 2.8. G <'^^ By for all y > 0. 

Proof. Indeed, since G F^o, G E'p, then, because of Fy<i <<^^ Fy, we also 
have G Fy for all y > 0. Therefore, 

(2.28) ryG + (1 - ry)G ryFy + (1 - ry)E,,. 

The left-hand side of (I2.28P is G'(x), while the right-hand side is By{x). Therefore, 

G By. □ 



Property 12.81 is crucial. Their consequence is the following result G Bn,k, 
where Bn^k{x) = Pi'{i?n,fc < a^}, which follows immediately from this Property 
by application of the formula for the total probability. Indeed, 



Bn,k{x) = Pr{i9„,fc < x} 
and consequently for j = 0, 1, . . ., 



/•OO 

/ By{x)d¥r{Tn.k < y}, 
Jo 



POO POO / POO \ 

J e-^"x^B„,fc(a:)dx = J e'^'^x^ i^J By{x)dPi-{TnM < y}) dx 

nOC p OO 

< / e-^"a;^G(x)da; / dPr{T„,fe < y] 
Jo Jo 

e~^'=x'G{x)dx. 



Hence, G Bn,k- 

Therefore, going back to (|2.23p we conclude that, according to Lemma Al of 
Appendix A, there exists the minimal random variable <r in the sense of the sto- 
chastic order relation, i.e. ^ <st Kn,k for all k = 1,2,..., corresponding to the 
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minimal probability distribution function G{x) in the sense of the order relation 
Specifically, 



(2.29) Pr{^ = m}^ /"°° e-^"^^^dG(x). 

Jo "i! 

Thus, taking the sequence of independent identically distributed random variables 
Ci, i = 1,2, all having the same distribution as i^, one can conclude 

/(") 

Therefore, the statement of Theorem 12.41 follows due to stochastic comparison. □ 



Corollary 2.9. In the case P — \ under Conditions (A) and (D) we have 

Zn. 



Zn denotes the number of offspring in the nth generation of the Galton-Watson 
branching process with Zq — I and the offspring generating function 

B{\) „ l-z + 2zB(X) , , , 

(7 z < ^—^ + 2e z <l 

1-z + zBiX) l-z + zB(X) 

where 



B{X) ^ pFiX) + {1 - p)- 

^ + A 

F{X) = / e^^^di^(x). 



?i, i = 1,2, . . is the sequence of nonnegative integer-valued independent identically 
distributed random variables all having the probability law 

Pr{Q - m} = / e-^^^^dG(x), 



and G{x) is a probability distribution function that satisfies the inequalities G 
FyO and G Efj_, where FyO <<^^ Fy for all y > 0. 



Proof. As we can see the only difference between the statements of Theorem 12.41 
and Corollarv l2.9l is in the coefhcient before e. Specifically, similarly to (|2.20p . 

/•OO 

fc(x) - B{x)\ < / sup \By{x) - dPr{Tj- fc < y } < 2e, 

Jo x>0,y>0 

V ' 

<2e due to f2M 

and then similarly to (|2.22p 

|B,,j(A) -i?(A)| < 2e 

as well, where 

=Pr{i?,-fc <x}. 

The rest part of the proof of the corollary is similar to the proof of Theorem l2.4l □ 
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3. The number of losses under Condition (B) 

In this section it is assumed, additionally to Condition (A) that the probability 
distribution function F{x) belongs either to the class NBU or to the class NWU. 
Having this additional assumption, Condition (B) enables us to obtain stronger 
inequalities than in the previous section under Condition (A). 

The lemma below gives a stronger result than characterization Lemma 12.21 
Specifically, we have the following lemma. 

Lemma 3.1. Under condition (B) we have 

,je{F, E^) < e. 

Proof. The proof of this lemma is given in Appendix B. 

Lemma 3.2. Under Condition (B) we have 

sup \By{x) — B{x)\ < 3e. 

x>Q,y>a 



□ 



Proof. We start the proof from representation (|2.2p . Then, as in the proof of Lemma 
2.11 we have (j2.3p . However, taking into account that F{x) belongs to one of the 
aforementioned classes NBU and NWU, then instead of the earlier inequality (|2.4p 
given in the proof of Lemma 12.11 we have the stronger inequality 

(3.1) \F{x)-e-n<^, 

which in turn is the result of the application of Lemma [3TT1 From (j2.2|) . (|2.3p . (|3.ip 
and the triangle inequality, for all x > 0, y > we obtain: 

\pF{x + y)~ p'F{x)F{y) - p{l - p)e-t'''F{y) \ 

F{x + y) 



< pF{y) 



< pF{y) 



_F{y) 

F{x 



+ p^F{y)\F{x) -e" 



(3.2) 



F{y) 
p''F{y)\F{x) 



-F{x)+F{x) 



-fiX 



< 



pF{y)[ 



\Fix) 



F{y) 

+ p'F{y)\F{x)~e-^'''\ 

< 2epF(y) + ep''F{y) 

< iepFiy), 

1(1 -p)e-^(^+^' -p(l -p)e-'^^F(x) - {I - pfc-^'^'^+y'^] 

< (1 -p)e-^«|e-'^^ - F(a;)| + {I - pfc-''«\c-'''' -F{x)\ 

< e(l - p)e-^'y + e(l - pfe-t'y 

< 3e{l~p)e-''y. 

Therefore, from (|3.2p and (|3.3p for all a: > 0, j/ > we have 



(3.3) 
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Lemma 13.21 is proved. 

Remark 3.3. In the special case where p ~ ^ wc have: 



□ 



(3.4) 



+ y)~ \F{x)F{y) - \e~^^-F{y)\ 



< 



\F{y) 
+ \F{y) 



< hF{y) 



F(y) ^ 
F{x+y) 



F{V) 



Fix] 



F{y)\F{x)-e- 



fix 



< hF{y), 



(3.5) 



lp-fj,ix+y) 



4^ 



F{x) 



4*- I 



Therefore, from (13.41) and (13.511 we have 



\By[x)-B{x)\ < ^ 



eFiy) 



F{y) 



Fix)\ 



He 



The main result of this section is the following 
Theorem 3.4. Under Conditions (B) and (D) we have 



Zn denotes the number of offspring in the nth generation of the Galton-Watson 
branching process with Zq ^ I and the offspring generating function 



9{z) < 



B{X) 



1- z + zB{\) 



3e- 



1 - z + 2zB{\) 
+ zB{X) ' 



\z\ < 1, 



1 



The sequence of independent identically dis- 
1,2,.. is as in Theorem\ 



where B{\) is as in Theorem 
tributed random variables <;i, i 

Remark 3.5. If Condition (D) is not satisfied, or the probability distribution Fyo{x) 
is unknown and cannot be evaluated, then a probability distribution function G{x) 
should be chosen such that G <<^^ Fy for all y > and G E^. 



Proof. The proof of the theorem repeats the proof of corresponding Theorem 12.41 
The only difference involves using the estimate given by Lemma 13.21 rather than 
estimate given by Lemma [2. II □ 

Corollary 3.6. In the case p — ^ under Conditions (B) and (D) we have 

z^ 
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Zn denotes the number of offspring in the nth generation of the Galton-Watson 
branching process with Zq ^ I and the offspring generating function 

B(X) 3e 1-Z + 2zB(X) , , , 

g{z) < + — ^r^, z < 1, 

l-z + zB{\) 2 l-z + zB{\) 

where B{X) is as in Theorem \2.4\ The sequence of independent identically dis- 
tributed random variables <;i, i = 1,2, is as in Theorem\2.4\ 



The proof of Corollary 13.61 is similar to the proof of Theorems 12.41 or 13.41 or 
Corollarv l2.9l The only difference is that the proof of this corollary uses the estimate 
obtained in Remark 13.31 rather than the estimate of Lemma 13.21 



4. The number of losses under Condition (C) 

In this section we study the number of losses under Condition (C). Applying Lemma 
13.11 we have the following. 

Lemma 4.1. Under Condition (C) 

je{B,E^) < e. 

The main result of the section is the following theorem. 

Theorem 4.2. Under Conditions (C) and (E) we have 



(4.1) ^^^<stLn<stY^ 



^ Vi. 

i=l 



Xn denotes the number of offspring in the nth generation of the Calton- Watson 
branching process with = 1 and the offspring generating function 

, , , , BiX) I- z + 2zB(X) , , 

4.2 9xiz)< + € 2 < 1, 

l~z + zBiX) l-z + zBiX) 

and Yn denotes the number of offspring in the nth generation of the Calton- Watson 
branching process with Yq = 1 and the offspring generating function 

gY{z) = , \z\<l. 

l-z + zB{X) 

The sequence of independent identically distributed random variables i = 1, 2, . . 
is determined as follows: 

(4.3) PrU,^m}= / e-^^-^^-dB„o(a;), m = 0, 1, . . . , 

Jo ml 

where ByO (x) is the minimal probability distribution function in the sense that 
ByO By for all y > 0. The sequence of independent identically distributed 

random variables Vi, i = \,2, . . ., in turn is determined as follows: 

(4.4) Pr{u, = m}= / e^^=^-^^dB(a;), m = 0, 1, . . . 

Jo ml 

Remark 4.3. If the probability distribution Byo(x) does not exist, then we only 
have the one-side stochastic inequality i„ <st X^i^i ^i- 
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Proof. Notice that under the assumption that B{x) belongs to the class NBU, 
we have the inequality J2i=i'^i —st Ln, where Vi, i = 1, 2, . . ., is a sequence of 
independent identically distributed random variables, defined by the probability 
law Pr{u, =m} = e-^^^dBix). 

Indeed, under this assumption from (|2.21[) and (|2.24p we obtain that Kj,k <st '«o,i 
for all < j < n (n > 0) and fc > 1 (for details see Abramov ,100, and then the 
desired inequality J2i=i ''^i —st Ln follows. 

In order to prove the second inequality L„ >st X]i"i '^i note the following. 

Assumption l\l.'S\i means that for all < j < n (n > 0), /c > 1 we have J^{'dj^k, x) < 
e. Then using the same arguments as in the proof of Theorem 12.41 we arrive at 
the left-side of inequality (|4.ip with the offspring generating function defined by 
inequality (|4.2p . The remaining part of the proof of the theorem is similar to the 
corresponding part of the proof of Theorem 12.41 □ 



Appendix A: The order relation 

Let ^1, ^2 be two nonnegative random variables, and let Si(a;) and S2(x) be their 
probability distribution functions respectively. 

Definition Al. The random variable is said to be less than the random variable 
^2 in the sense of the relation (notation: <<^^ ^2 or Si <<^^ S2) if for fixed 
parameter A > and all i = 0, 1, ... there is the inequality 



/■oo /"C 

/ e-^^x'Si(x)dx > / 
Jo Jo 



x'^'E,2{x)dx. 



Relation is a partial order relation. It follows from the following lemma. 

Lemma Al. Let ^1 <'^^ ^2- Then for the random variables 61 and 62 given by the 
probability laws: 

Pr{0fe = m} = / c-^-^^dSfc(a;), A: = 1, 2; m = 0, 1, 
Jo rn\ 

we have 9i <st 82- Converse, if Oi <gt 62, then ^1 <<g'^ ^2- 
Proof. Indeed, by partial integration we have 

POO /"OO 

/ e-^MSfc(x) = A/ e~^"Sfc(a;)dx (fc = l,2), 
Jo Jo 



and 



-d^fc(x) 



TO 

= A / e — '—^k{x)dx~ / e — '——:ik{x)dx 
\Jq to! Jo (m~l)! 

(to = 0,1,...; k = 1,2). 

Therefore, for all to = 0, 1, and fc = 1, 2, we obtain 

A re-'^^-^E,ix)dx^f2 re~^^^dS,(x), 
Jo ml j^Jo tl 

and the result follows. □ 
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The only difference between the proofs of these lemmas is that there is an addi- 
tional condition (B) in Lcmma l3TI Therefore, we provide the general proof of both 
Lemmas 12.21 and 13.11 with an appropriate specification of the case studies. 

We have 

F{x + y)^F{v)F{x) + h{x,v), 
where y) = F{x) — Vv{C, < x + y \ Q > y}. Then, 



{B.l) 
where 
(B.2) 
Denoting 



1 — 

F{x + y)Ay = -F{x) 



1 



e(x). 



e(x) = M/ F{y)h{x,y)Ay. 



^{x) 



F{x + y)dy^ / F{y)dy, 



from (B.l) we obtain the following differential equation 
(B.3) *'(a;) - e(x) = 0. 

The solution of differential equation (B.3), satisfying the initial condition ^'(0) 
I/H, is 



1 



Hence, 
(BA) 



F{x) =e-^^+^ / e 



Q{t)dt. 



Q{t)dt~Q{x). 



Let us study some properties of (B.4). Note first that |0(a;)| < e. Indeed, according 
to assumption of this lemma sup3.>o^j^>o \h{x,y)\ < e, and therefore, as it follows 
from (B.2), 



le 



{x)\ < ( sup \h{x,y)\ ]fi 

^ x>Q,y>0 



Therefore, 



{B.5) 



-At(a;-t)p 



Q{t)dt-Q{x) 



< 






Jo 



F{y)dy < e. 



< 2 sup \e{t)\ < 2e. 

0<t<x 



Substituting (B.5) for (B.3), for all a; > we obtain 

|F(x) -e-'^^l < 2e, 

and the statement of Lemma [2.21 follows . 

Under the assumption that F{x) belongs either to the class NBU or to the class 
NWU, the function Q{x) is either nonnegative or nonpositive. This is because 
h{x, y) contains the term belonging to one of these classes, and therefore is either 
nonnegative or nonpositive. Therefore, according to (B.2) 



{B.6) 



-t-i^~t)e{t)dt - e{x) < \e{x)\ < e. 
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Substituting (B.6) for (B.3), for all a; > we obtain 

\F{x)-e-^^\ < e, 
and Lemma l3. II is therefore proved as well. 

Appendix C: The proof of Lemma 

Let ^ := {tn,n > 0} denote the time instants at which customers arrive to 
the system (a customer can arrive without joining queue), where tg = 0, and let 
{N{C),C S ^{M.)} denote the Poisson process that consists of the collection of 
points i^. 

Next, let {T„, n > 1} denote the time instants at which a new service begins in 
the system. At these time instants either 

(i) a departure occurs and leaves at least one customer behind the system, or 

(ii) an arrival occurs and finds an empty system. 

Finally, let Q{t) denote the number of customers in the system at time t. 
Let us find Pr{i9j_fe < x\Tj^k = u}- Suppose ^„ = <Ar(T„)+j-Q(T„) - T^- Clearly, 
is the remaining amount of time, after the beginning of the mth service, until a 
new customer arrives, plus the next j — Q{Tm) — 1 interarrival times. The instant 
Tm is chosen such that Sj,k-i < Tin < ^j,fe, where we set sj^q := 0. (The definition 
of the time instants sj^h and tj,k is given in the proof of Theorem 12.41 ) Denoting 
r] = inf{n > 1 : Xm > £.7n}, where Xm denote the length of the mth service time, 
one can see that 

Pr{i?i,fc < x\Tj^k = J/} = Pr{x^ - C»? < 2^1^^ = y}, 

and 

Pr{x-n -£,v < x\U = y} 

oo 

= ^ Pr{x,, - Crj < x\t] = n, f„ = y}Pr{77 = n\C^ = y} 

n=l 

oc 

= ^ Pr{xn -^n< AXra < £,ni,l < m < n - 1, Xn > V , ^ y}Pr{r] = n|^^ = y} 

71=1 

^Pr{Xn<X + y,Xn>y\Xm<£,rnA<m<n-l,^ri=y}^ f . 

Notice, that each C„ is a function of the Poisson arrival process and the random 
variables xi, X2, • • ■ , Xn-i, and therefore the random variable %„ is independent 
of the vector (^i, xi, 6, X2,- • ■ , C„), and so 

Pr{Xn <x + y,Xn> y\Xrn < £,m-A < m < Tl ~ 1 , C„ = y} 

Pr{Xn > y\Xm < Cm, 1 < "1 < 71 - 1,C„ = y} 

Hence, 

Pr{x>7 - C»? < x\^v =y} = By{x), 

or 

Pr{-dj^k < x\Tj^k = y} = By{x). 
and the statement of Lemma 12.61 is proved. 



Byix). 
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